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ABSTRACT 


A fast method of retrieving vertioal temperature profiles in the 
atmosphere and of determining the paths of the rays producing the ATMOS 
oocultation speotra has been developed. 

The results from one set of oocultation data appear to be consistent with 
other available data, 

A study of sources of error, a search for other suitable features for 
measurement in the speotra, and modification of the program to obtain mixing 
ratio profiles have been initiated. 
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INTRODUCTION 


One of the principal objectives of the ATMOS experiment is to obtain 
vertical mixing-ratio profiles (VMRs) of atmospheric gases. This requires the 
geometry of the ray paths and the vertical temperature and density profiles in 
the atmospheric region sampled. This information can be obtained from sources 
external to the experiment but:, its quality may not be sufficient to allow the 
highest accuracy in the VMR profiles of interest. 

This information can be retrieved from the occultation spectra and wc 
have demonstrated that temperature profiles and ray paths can be obtained from 
the ATMOS data quickly and precisely. However, in order to analyze all the 
available spectra, additional absorbing features must be identified and a 
thorough error analysis completed. 

Descriptions of the method, the results and the current stage of error 


analysis are given. 



METHOD 


The analytical method has been described by Hoke and Shaw (Ref. 1). A 
series of suitable lines of gares with known vertioal mixing ratios and with 
accurately determined line parameters are identified in the occultation 
spectra and their equivalent widths measured. In our retrieval of the 
temperature profile from sunset occultation number 6, obtained on April 30, 
1985, up to 30 lines in sixteeen spectra when chosen for analysis. These 
provided more than 120 equivalent widths. Some of these lines originated from 
low lying states and had temperature insensitive intensities. The equivalent 
widths of these lines allow the total amount of absorber along the path and 
hence in the ray path itself to be estimated (1). Another series of lines 
arising from high energy states was chosen beoause these lines intensities are 
strongly temperature dependent. 

A nonlinear, least-squares (NLLS) program was then used to adjust the 
initial guesses for the temperature profile and the tangent heights of the 
rays until a best fit between the calculated and observed equivalent widths 
was obtained. 

The parameters retrieved included the temperatures at fixed altitudes In 
the atmosphere and five parameters describing the dependence of the tangent 
heights of the rays on time. 

In order to obtain a retrieval it is necessary to make numerous 
assumptions. These include: 

1. The line parameters (intensity, position, shape, width, and lower state 
energy are accurately known. 

2. The atmosphere is in hydrostatic equilibrium. 

3. The VMR of C0 2 is known. 
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4. There are no systematic errors in the equivalent widths. 

5. The effects of winds can be included. 

6. The effeots of LTE/non LTE can be modelled. 

7. Various smearing effeots caused by the finite FOV of the instrument, by 
satellite motion, and by jitter in the sun tracker oan be inoluded in the 
analysis. 

8. The atmospherio pressure at some reference altitude is known. 

The effeots of errors in some of these assumptions is discussed below. 
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RESULTS 


The temperature profile in Fig. 1 was obtained by analyzing over one 
hundred equivalent widths of thirty lines of C0 2 ooourring in sixteen spectra 
obtained during a sunset ocoultation observed by ATMOS on April 30, 1985. 

The tangent heights of the rays lay between 15 and 85 km. Eleven 
temperatures at altitudes between 20 and 80 km were retrieved and these are 
compared with a temperature profile supplied by Oxford University in Fig. 1. 
By analyzing subsets of the data and by choosing different spacings and 
altitudes for the temperature parameters it has been shown that the observed 
features are reproducible. 

Less than three minutes of computer time were required to obtain these 
results. 
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The temperature profiles retrieved by using different heights for 
the temperature parameters. 





SOUnCES OF ERROR 


Line Parameters 

The lower 3tate energies of all the C0 2 lines used in this work are 
aoourately known and the use of equivalent widths does not require an aoourate 
knowledge of the line positions, The doppler widths e D oan be oalo.ti.ated if 
the temperature is known. There are uncertainties in the line intensity S, 
the line shape, and the oollision-broadened widths a L . 

Intensity Errors 

The equivalent widths of the C0 2 lines analyzed do not lie in the linear 

region of the curve of growth. Sinoe the linns lave essentially a Voigt shape 

the variation of the equivalent widths on absorber amount is a function of 
a D /aL and hence of the height of the ray in the atmosphere, We have found 
that fractional errors AS/S of 5 % in the temperature insensitive lines chosen 

cause errors in the estimated tangent heights of the rays of 0.^ km -»r less at 

all altitudes. 

Fractional errors AS/S of 535 in the intensities of the temperature 
sensitive lines oause errors in the temperature profile which depend on E" and 
on the position of the equivalent width on the curve of growth. For the lines 
chosen errors of less than IK occur for lines with E" > 3000 cm - ' but this 
increases to more than 3K for lines with E M £ 800 cm -1 . 

Line Widths and Shapes 

Since a D /a L << 1 for lines formed above 50 km and sinoe the lines chosen 
for measurement do not have significant contributions to the equivalent widths 


from the Lorentz wines It is not expected that uncertainties in the collision 
broadened widths will significantly affeot the results. It was found that 
increasing the width a L by 2055 ahanged tne retrieved temperatures by less than 
2K between 20 and 80 km altitude. 

Volume Mixing Ratio 

The C0 2 volume mixing ratio is well known below 60 km. At higher 
altitudes dissociation ocours and there i3 considerable uncertainty in this 
VMR profile. 

We have found that a change of 255 in the C0 2 VMR causes changes of les3 
than 0.2 km in the retrieved tangent heights and the change in the temperature 
profile is less than IK. Muoh larger uncertainties in the VMR above 60 km are 
anticipated. In this oase it may be neoossary to assume the ray paths are 
known and to retrieve the temperature and VMR (C0 2 ) simultaneously. We have 
not yet explored this problem in depth. 

Reference Pressure 

In order to describe the ray paths through the atmosphere, obtained from 
the spectra, in terms of the geometry of the oooultation, it is necessary to 
assume that the atmospheric pressure at 3ome referenoe altitude is known. 

This problem has been discussed by Shaffer et al (2) who have shown that a 
fractional error AP/P of 5/5 in the referenoe pressure causes a corresponding 
change in the tangent heights of the rays of less than 0.5 km. 

However, errors in the tangent heights of the rays introduce errors in 
the temperature profile. We have found that an error in the tangent height of 
1 km, which corresponds to an error in the reference pressure of 10^, causes 
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an error In the temperature parameters which increases with height above the 
reforenoe pressure level from about 3K/km near 30 km to 6-7 K/km near 80 km. 

Winds 

Variations in the component of the wind along the line of sight and at 
different levels in the atmosphere oau3e line shifts, These shifts are small 
but should be observable in the ATMOS 3peotra <see Appendix A), Provided the 
lines are in the linear region of the ourve of growth these wind speed ohanges 
will not alter the equivalent widths of the lines. However, as the lines 
saturate, wind speed differenoes cause the observed equivalent widths to 
increase. We have found that, provided the wind speed is less than 100 m/s 
the maximum fractional increase in the equivalent width iW/W is less than 3? 
at all heights and at all frequencies, This is only slightly larger than the 
estimated uncertainty in measuring the equivalent widths of the lines and is 
therefore not eona.iri.r-’od to be an important effect. 

Spectral Smearing 

Because of satellite motion during the period of data collection, the 
finite field of view (FOV) of the instrument, and Jitter in the sunseeker the 
collected 3peotra are some weighted mean of many atmospheric paths. For a 2mr 
FOV these effects are equivalent to a spread in the tangent heights of the 
rays reaching the instrument of more than 6 kin. If speotra are averaged this 
spread is increased. 

These effects limit the vertical resolution of the temperature and VMR 
profiles and they may also affeot the accuracy of the retrieved profiles 
unless they are included in the data analysis. In particular, if there is a 
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significant vertical temperature gradient over the FOV then the absorber-mas 3 - 
weighted mean tangent height determined from the pressure lines may differ 
from the line-intensity-mass-weighted mean tangent height of the temperature 
sensitive lines. Then effeots have not yet been thoroughly explored. It Is 
expeoted that soma method of weighting for these effeots must eventually be 
Included in the calculations. One method is evaluate mean speotra by 
considering the contributions of several ray paths to each spectrum. 

Non-LTB Effeots 

Departures from LTE are expected to ooour in the atmosphere above about 
50 km. We have not yet Included these effeots in our calculations or con- 
sidered uncertainties in these effects on the retrieved temperature profiles. 

Sensitivity 

In the previous section we have discussed the effeots of systematic 
errors in the assumed models and model parameters on the retrieved temperature 
profiles. We have found, by changing these models and model parameters, that 
significant changes in the retrieved vertical temperature profiles can be 
made. 

However, in all the cases we have explored, it has been found that the 
solutions are strongly convergent and that the final profiles are obtained 
after only a few interations. 

Although, we expect significant effects due to improvements in the COg 
VMR and to the inclusion of non-LTE effects in the temperature profiles 
obtained at 100 km and above it appears that the precision with which most of 
the other required parameters is already known is sufficient to obtain 
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temperature profiles accurate to 2-3K and’hence that the principal souroe of 
uncertainty is due to the information content of the spectra. 

This information is limited by the SNR of the collected spectra, by the 
frequency spaolng of the speotral data, by the finite iistrumental resolution, 
and by the necessity for using partially saturated lines, since the equivalent 
widths of lines in the linear region of the curve of growth are too small to 
be measured accurately. 

For lines in the linear region of the curve of growth the ratio (AW/W)/ 
(Au/u) when u is the absorber amount is unity. The corresponding value of 
thl3 ratio for the COg lines measured in this study was, in some cases, as low 
as 0,2. 

Despite this laok of sensitivity and the measurement precision of about 
2 % in the equivalent widths the precision of a few degrees in the retrieved 
temperature profile is excellent at this stage in the analysis. 
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Summary of systematic Error Effeots 


These estimates are based on the analysis of simulated data similar to 
those obtained from ATMOS oooultation SS06. 


Effeot 

Magnitude 
of error 


Effect 
on h T 

Effect on 
T profile 

Measurement of 
equivalent width 

AW/W - 5 % 

< 

0.7 km 

1 .5 - 5K 

Line Intensity 

AS/S - 5 % 

< 

0.5 km 

1 - 3K 

(depends on E") 

Collisional width 

Act/cx ■ 20? 

< 

0.5 km 

< 2K 

Volume mixing ratio 

A(VMR)/VMFl - 2 % 

< 

0.2 km 

< IK 

Pressure at 
reference altitude 

AP/P » 5? 

< 

0.5 km 

3 - 8K 

Winds 

AW/W < 3? 
for v <* 100 m/s 

< 

0. 4 km 

1 - 3K 

Spectral smearing 


< 

B 

in 

o 

1.5 - 

LTE/Non-LTE 


to be investigated 
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Abstract 

By measuring the mean shifts of hundreds of lines in solar occultation 
spectra obtained at satellite altitudes relative wind speeds along the lines 
of sight of the rays can be obtained with precisions of S m/s or better at 
altitudes up to 100 ton. It is necessary to obtain more accurate line 
positions or to Introduce a gas sample for frequency calibration if absolute 
wind speeds are to be obtained. 
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Introduction 

A series of solar occultation spectra is expected to be obtained with a 
Fourier Transform Spectrometer** carried on the Space Shuttle, This ATMOS 
instrument collects spectra in broad regions between 2 and 16 im> at a rate of 
about one per second and with a spectral resolution of about 0.01 cm" 1 as the 
tangent heights of the rays vary from 0 to over 150 km. Solar lines occur in 
spectra of large tangent heights and ( as the rays penetrate more deeply into 
the atmosphere , thousands of lines belonging to more than thirty atmospheric 
gases are also formed. Many of these lines are due to COg which has strong 

bands near 2.0 t 2.7, 4.3 and 15 pm. Some lines of COg appearing in a simulated 

spectrum with a resolution of 0.01 cm" 1 and a tangent height of 85 km are 
shown in Fig. 1. 

Because of the relative speeds v of the satellite with respect to the 
undisturbed solar atmosphere or the earth's atmosphere along the lines of sight 
of the rays the spectral features are shifted by amounts Av where 

Av/v - v/c. (1) 

These shifts can be calculated from the orbital geometry parameters. 

Atmospheric motions or wind components along the line of sight Introduce 
additional line shifts. By measuring these shifts as the FOV moves over the 
solar disk information about motions in the solar atmosphere can be obtained. 

Winds of up to 100 m/s occur in the earth's atmosphere. These speeds vary 

with altitude and location and hence telluric lines with contributions from 
different levels in the atmosphere may have distorted shapes and shifts. 

“The ATMOS instrument is primarily intended for investigating the composition 
of the stratosphere and mesosphere. The principal investigator is C. B. Farmer 
of the Jet Propulsion Laboratory. 
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Useful meteorological information can bo Inferred if the wind speed components 
along the line of sight as low as 5 m/oec (10 oph) can be measured. These 
components are described as winds in the remainder of this discussion where it 
is shown that these precisions can be achieved if the Instrument is stable. 

The wind information obtained from the ATMOS spectra has a high vertical 
resolution over the region from below 30 km to above 90 km. A number of other 
methods for determining winds in this region by satellite remote sensing have 
been proposed. All are based on detecting the effects due to Doppler shifts 
and they have been summarized or referenced in papers by Reese (1), Hays (2), 
and McCleese and Margolis (3). Some techniques (1,2) measure the shift of a 
single line observed with a high resolution Fabry-Perot interferometer or the 
combined effects due to many lines as measured with a ncn-dispersive 
instrument (3). The precision of the method discussed below is enhanced by 
the measurement of the mean shifts of hundreds of lines. 

Line Shapes and Shifts 

Most absorption lines in the earth's atmosphere can be described by the 
Voigt shape. In the troposphere this approaches the Lorentz shape with 
typical widths of somewhat less than 0.1 cm~ 1 /atm. Above the tropopause the 
shape approaches the Doppler shape and the widths are given by 

c^/v - 7.16 x 10" 7 (T,'M)*, 

where T is the temperature and M 1 b the molecular weight of the absorber. 

Over the range 0 to 100 km atmospheric temperatures vary from about 200 
to 300 K and the molecular weights of many absorbing molecules lie in the 
range from about M - 16 (CH^) to M - 44 (NgO, CCg). Doppler widths for 
several seta of values of T and M over the region from 500 to 5000 cm -1 are 
given in Table 1. 
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The Doppler shifts caused by winds of 5 m/a and 100 m/a over this same 

region are shown in Table XI. Those shifts are small compared with the line 

widths in Table 1 and also with the resolution K of about 0.015 cm" 1 of the 

apodized (4) ATMOS spectra. Hawkins ot al (5) have found that the positions 

of lines with widths narrower than the spectral resolution, observed in Fourier 

transform spectra with SNR's in excess of 100 can be estimated with precisions 

of R/30 or better. The ratios P • 30Av<5 a/s)/R in Table II show that, with 

this criterion, the minimum wind speed of 5 m/a cannot be estimated from a 

single line. However, the required precision can be achieved by measuring 
2 

N - 1/P line positions (5). The value of N in Table II increase rapidly as 
the line frequency decreases because the Doppler shift is directly proportional 
to frequency. 

The shifts caused by a 5 a/s wind are smaller than the accuracies with 
which the frequencies of many lines of gases such as C0g, H^O, N 2 0, CO, or CH^ 
are known <3). Thus, although it is possible to measure the relative shifts of 
groups of lines due to winds in different spectra, and hence to obtain the 
relative wind speeds, the determination of the absolute wind speed from the 
spectra requires additional spectroscopic measurements. Until these are 
available it will be necessary to normalize the relative wind speeds by direct 
comparison with in situ wind measurements. 

Choice of Lines 

The lines occurring in the spectra include: 

1. Solar lines. These appear in all the spectra but many are obscured 
by telluric lines at low tangent heights. These solar lines do not yield 
direct information about motions in the earth's atmosphere. However, by 
measuring the relative shifts of groups of lines observed in successive 
spectra of large tangent heights, the short term stability of the 


Instrument and the guiding accuracy of the sunseeker Can be inferred. The 
changes in the relative speed of the satellite and the solar atmosphere 
over the portion of the solar disk viewed can also be determined. 

2. Lines of stratified atmospheric gases. Some gases such as 0^ have 
strongly peaked mixing-ratio profiles and produce largo numbers of 
measurable lines in the spectra. Below some critical height the peak 
absorptances of these lines increase until the layer of maximum concen- 
tration is penetrated. If this Increase is sufficiently rapid the wind- 
induced shifts of the lines in each spectrum are dominated by the windB 
in the lowest layer penetrated. Also, if there is little absorber in 
the atmosphere below the concentration maximum, the peak absorptances of 
the lines in spectra obtained for lower tangent heights no longer 
Increase and the wind shifts are frozen. Such lines can be used as 
markers with which new lines of other gases appearing in spectra below 
this peak can be compared. 

3. Lines of gases with nearly constant mixing ratio profiles. Although 
the peak absorptances of lines of gases such as C0 2 depend on the line 
Intensity and spectral resolution they all Increase steadily as the 
tangent height decreases. An example of this growth for a line of a gas 
with a constant mixing-ratio in the standard atmosphere is shown In 
Pig. 2. The dependences of the peak absorptances of other lines with 
different intensities on tangent height are shown in Pig. 3. Since line 
positions can usually be most precisely measured if the peak absorptances 
lie In the range 20 to 80% these curves indicate that each line can be 
used over a range of about 20 km and that several sets of ’ lines are 
required to cover the region of interest. 
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Estimation of Wind Speeds from Lines of a Gas with a Constant Vertical 
Mixing-Ratio 

The relative velocity of the satellite and the "wind free" atmosphere is 
essentially the same for all the molecules along the atmospheric path and thus 
it produces the same shift for all the "components" of the line formed in the 
30 km deep layer where most of the absorption occurs. Uncertainties in tlv» 
satellite height of 10 km or of 1° in the latitude of the closest approach of 
the ray to the earth cause errors of less than 5 m/s in this motion. Since 
the orbital parameters are known to this accuracy or better the shift due to 
satellite motion can be accurately calculated and the absolute wind speeds can 
be obtained from the residual shift, provided the absolute frequencies of the 
lines are known. 

If the entire atmosphere above some level h^ has a wind speed' v^ then 
all the components of a line in spectrum S, observed for rays with a tangent 
height hj^ are wind shifted by 

*»1 - *1 - *0 • 

” v i v c • <a| 


where v is the position of t&£ line in the absence of a wind and A\>, « v . 
o l o 

The contributions by the absorber above h^ to the same line in spectra with 
lower tangent heights will also be shifted by A\>^. In spectrum Sg, obtained 
for a ray path with a tangent height hg " h l " &h 2 * fl ame line has an 
increased peak absorptance due to the absorber AUg along the path in the layer 
Ahg. If there is a wind Vg in this layer the contribution to the line by AUg 
in spectrum Sg and other spectra with tangent heights leas than hg is shifted by 

4 ”a ■ *2 " *o • 


- v 2 V° ' 
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Above 30 km Unas which are not strongly saturated have shapes which are 
nearly invariant with height and the not shift Av 2 of the observed line in 
spectrum S 2 with respect to con be considered, to a first approximation, 
as the mean of the line shifts weighted by the amounts of absorber above and 


below h 1 
In general 


4v 2 « (Au^ + AUgVgJv^UUj^ + Au 2 )c. 


♦* 

n n 

*v n « < v Q /c) l AUjVj/ J Au, , 

i -1 1-1 


If, for example, all the winds are zero except in the bottom layer Ah R the 
shift 

l» n - <V c > 4 VnV u • 

where 

n 

u - l Au A , 

1»1 

Hence 

4v n /Av n * Au n /u * (3) 

As the thickness Ah increases AU increases and Av Av . 

n n n n 

Although the line shift can be measured more precisely the vertical resolution 
decreases. 

When Au n « u/2 the ratio of the shifts in Kq. 3 is a. 5. the change in 
tangent height required to double the amount of an absorber, with a constant 
vertical mixing-ratio, along the ray path varies from about 4.5 to 5.5 km over 
the altitude range 0 to 100 km. Thus, if it is desirod to obtain the winds in 

layers 5 km thick by measuring the shifts of COg lines, these shifts are 

/ 

approximately one half those given by Bq. (2) and four times the number of 
lines given in Table II must be measured to obtain a precision of 5 m/a. 


8 


Inspection of lists of the parameters of atmospheric lines and of simulated 
solar spectra such as that In Fig. 1 Indicates that there are sufficient 
numbers of strong lines near 2.7 and 4.3 urn to allow winds to be measured up 
to 100 km. 

Discussion 

It has been shown that winds speeds can be determined with accuracies of 
5 m/s from ATMOS spectra provided the orbital parameters, and the line 
positions arc accurately known and the Instrumental performance is 
reproducible. Ground based tests have Indicated that the Instrument Is stable 
and the stability in orbit can be confirmed by measuring the positions of 
solar lines. These solar lines can be used to transfer the frequency 
calibration from one spectrum to another provided that motions in the solar 
atmosphere are small. In future experiments the system con be internally 
calibrated with the method described by Hawkins et al (3) by introducing an 
absorption cell containing a mixture of absorbing gases. 

The number of lines to be measured to obtain the desired precision la 
only a small fraction of the total number of observable lines. Indeed the 
reproducibility of the results can be compared by measuring the mean shifts of 
subsets of lines selected on the basis of their peak absorptance, origin, 
separation from other lines, etc. 

In this discussion it has been assumed that the line shape is invariant 
with height and that the mean line shift la the mass weighted shifts of the 
individual layers. Better estimates of £he relations between the observed 
mean line shifts and the wind speeds can be obtained by computer simulations. 
Synthetic spectra covering wide spectral regions and incorporating the effects 
of winds, spectral resolution, finite FOV, and satellite motion can be 
calculated if the tangent heights of the rays, the molecular line parameters 
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Including the positions, intensities and shapes and the atmospheric character- 
istics such as the vertical pressure, temperature and absorber mixing-ratio 
profiles are known (7). If necessary, the effects of pressure shifts of the 
lines can be Included although HcCleese and Hargolls (3) note that this effect 
should be negligible above 40 las. 

Even though the vertical resolution is high the amount of Information 

which can be obtained from occultation spectra is limited by the low spatial 

resolution caused by the long atmospheric paths and because each of the 

% 

recorded spectra samples a different portion of the atmosphere. Thus only 
mean values of the atmospheric properties such as temperature, absorbing 
mixing-ratio or wind spoed can be inferred. 

In addition to affecting the line positions wind shears may also affect 
the shapes and the equivalent widths of the lines. Since most lines which are 
observed in spectra with tangent heights above 3G ka and which have abaorp- 
tances of 20% or greater lie outside the linear region of the curve of growth, 
the differential line shifts will Increase the observed equivalent widths. 

The effects caused by most of the winds encountered in the atmosphere are not 
likely to introduce significant orrors in, for example, the determination of 
the vertical mixing-ratios of trace gases. These effects soy, however, limit 
the accuracies with which isotopic abundances can be determined and also the 
pressure or temperature profiles derived from the spectra. Most of the methods 
for Inferring these profiles are based on the differences in the behavior of 
lines with temperature sensitive and insensitive intensities and are strongly 
dependent on the ability to fit the observed absorption accurately. Under some 
circumstances it may bo necessary to retrieve the wind field simultaneously 
with the temperature profile. 
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Tabic I 

Range of Doppler widths of lines of 
atmospheric molecules between 0 and 100 ka 
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Table II. Wind induced Doppler shifts , and the 
numbers of line positions required to measure 5 m/s winds 


Wave- 

number 

Doppler Shift fiv(cm -1 ) 

Resolution 

P - 

30Av(5 m/a) 

No, of lines 

(cm -1 ) 

V-5 O/S 

v-100 m/s 

R<cm -1 ) 

R 

N - 1/P 2 

500 

2X10 -6 

1.7X10 -4 

0.015 

0.016 

4000 

2000 

3X10 -9 

6.2X10 -4 

0,015 

0.06 

280 

5000 

8xl0 -9 

1.7X10 -3 



0.015 

0.16 

40 
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figure Captions 

Fig, 1 Absorption by COg In a simulated occultatlon spectrum with a 
resolution of 0.01 ca -1 and a tangent height of 85 km. 

Fig. 2 Simulated spectra of a line of a gas with a constant vertical mlxlng- 
ratlo observed with a resolution of 0.015 cm' 3, and tangent heights of 
43, 47, 51, and 55 km„ 

Fig. 3 The dependence of the peak: absorptance of lines of a gas with a 

constant vortical alxlng~ratlo on tangent height. These values were 
obtained from simulated spectra with a resolution of 0.015 cm -1 and 
the Intensities of the lines change by a factor of 10 In successive 
curves. The topmost curve corresponds to one of the stronger lines of 
COg near 15 pm. 
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